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ABSTRACT 

We investigate the environmental impact of the first active galactic nuclei that may 
have formed ~ 150 Myr after the big bang in low-mass ^ 10^ minihaloes. Using 
ENZO, an adaptive-mesh refinement cosmological hydrodynamics code, we carry out 
three-dimensional simulations of the radiative feedback from 'miniquasars' powered by 
intermediate-mass black holes. We follow the non-equilibrium multispecies chemistry 
of primordial gas in the presence of a point source of X-ray radiation, which starts 
shining in a rare high-a peak at z = 21 and emits a power-law spectrum in the 0.2-10 
keV range. We find that, after one Salpeter time-scale, the miniquasar has heated up 
the simulation box to a volume-averaged temperature of 2800 K. The mean electron 
and H2 fractions are now 0.03 and 4 x 10~^: the latter is 20 times larger than the 
primordial value, and will delay the buildup of a uniform UV photodissociating back- 
ground. The net effect of the X-rays is to reduce gas clumping in the IGM by as much 
as a factor of 3. While the suppression of baryonic infall and the photoevaporation 
of some halo gas lower the gas mass fraction at overdensities 6 in the range 20-2000, 
enhanced molecular cooling increases the amount of dense material at 5 > 2000. In 
many haloes within the proximity of our miniquasar the H2-boosting effect of X-rays 
is too weak to overcome heating, and the cold and dense gas mass actually decreases. 
We find little evidence for an entropy floor in gas at intermediate densities preventing 
gas contraction and H2 formation: we show, instead, that molecular cooling can affect 
the dynamics of baryonic material before it has fallen into the potential well of dark 
matter haloes and virialized. Overall, the radiative feedback from X-rays enhances gas 
cooling in lower-cr peaks that are far away from the initial site of star formation, thus 
decreasing the clustering bias of the early pregalactic population, but does not appear 
to dramatically reverse or promote the collapse of pregalactic clouds as a whole. 

Key words: cosmology: theory - galaxies: active ~ methods: numerical - X-rays: 
general 



1 INTRODUCTION 

The first 'seed' black holes (BHs) that later grew to become 
the supermassive variety that power active galactic nuclei 
(AGNs) must have appeared at very early epochs, z > 10: 
this is in order to have had sufficient time to build up via 
gas accretion a mass of several x 10^ M0 by 2: = 6.4, the 
redshift of the most distant quasars discovered in the Sloan 
Digital Sky Survey (Fan et al. 2003). The origin and nature 
of this seed population remain uncertain. Numerical simula- 
tions performed in the context of hierarchical structure for- 
mation theories show that the first stars (the so-called 'Pop- 
ulation Iir) in the Universe formed out of metal-free gas in 
dark matter 'minihaloes' of mass above a few xlO^ M0 
(Abel, Bryan, & Norman 2000; Fuller & Couchman 2000; 
Yoshida et al. 2003; Reed et al. 2005) condensing from the 
rare high-cr peaks of the primordial density fluctuation field 
at 2 > 20, and were likely very massive (e.g. Abel, Bryan, & 



Norman 2002; Bromm, Coppi, & Larson 2002; see Bromm & 
Larson 2004 for a recent review) . Non-rotating very massive 
stars in the mass window 150 < m* < 250 M0 are expected 
to disappear as pair-instability supernovae (Bond, Arnett, 
& Carr 1984) and leave no compact remnants. Stars with 
40 < m, < 150 Mq and m, > 250 M0 are predicted instead 
to collapse to BHs with masses comparable to those of their 
progenitors (Fryer, Woosley, & Heger 2001). Barring any fine 
tuning of the initial mass function of Population III stars, 
intermediate-mass BHs - with masses above the 5-20 Mq 
range of known 'stellar-mass' holes - may then be the in- 
evitable endproduct of the first episodes of pregalactic star 
formation (Madau & Rees 2001). Another route for the cre- 
ation of more massive black hole seeds may be the formation 
- and subsequent collapse as a result of the post-Newtonian 
instability - of supermassive stars with m, S> 10'^ M© (see, 
e.g., Shapiro 2004a for a recent review) out of the lowest 



© 2005 RAS 



2 Kuhlen & Madau 



angular momentum gas in rare haloes above the minimum 
meiss for atomic cooling at high redshifts (Bromm & Loeb 
2003; Koushiappas, Bullock, & Dekel 2004). Our simulations 
are related more closely to the former type of scenario. 

Physical conditions in the central potential wells of gas- 
rich protogalaxies may have been propitious for BH accre- 
tion. In the absence of an H2 photodissociating UV flux and 
of ionizing X-ray radiation, three-dimensional simulations 
of early structure formation show that the fraction of cold, 
dense gas available for accretion on to seed holes or star for- 
mation exceeds 20 per cent for haloes more massive than 
lO'' Mq (Machacek, Bryan, & Abel 2003). There is also the 
possibility that the collapse of a rotating very massive star 
may lead to the formation of a massive BH and a 'ready- 
made' accretion disc (Fryer et al. 2001; Shapiro 2004b) that 
may provide a convenient source of fuel to power a 'mini- 
quasar'. 

The presence of accreting BHs powering Eddington- 
limited miniquasars at such a crucial formative stage in 
the evolution of the Universe presents a challenge to mod- 
els of the epoch of first light and of the thermal and ion- 
ization early history of the intergalactic medium (IGM), 
and serves as the main motivation of this paper. Feed- 
back processes from the first stars and their remnants likely 
played a key role in reheating and structuring the IGM 
and in regulating gas cooling and star formation in pre- 
galactic objects. Energetic photons from miniquasars may 
make the low-density IGM warm and weakly ionized prior 
to the epoch of reionization breakthrough (Madau et al. 
2004; Ricotti, Ostriker, & Gnedin 2005). X-ray radiation 
may boost the free-electron fraction and catalyze the for- 
mation of H2 molecules in dense regions, counteracting their 
destruction by UV Lyman- Werner radiation (Haiman, Abel, 
& Rees 2000; Machacek et al. 2003). Or it may furnish an 
entropy fioor to the entire IGM, preventing gas contraction 
and therefore impeding rather than enhancing H2 forma- 
tion (Oh & Haiman 2003). Photoionization heating may 
evaporate back into the IGM some of the gas already in- 
corporated into haloes with virial temperatures below a few 
thousand kelvins (Barkana & Loeb 1999; Haiman, Abel, & 
Madau 2001; Shapiro, IHev, & Raga 2004). The detailed con- 
sequences of all these effects is poorly understood. 

In this paper, we describe the results of fully 3D Eule- 
rian cosmological hydrodynamical simulations of the effect 
of the first miniquasars on the thermal properties of the 
high-redshift IGM. The focus of our investigation is not the 
modeling of the processes that lead to gas accretion on to 
BHs, but rather the radiative feedback of this population on 
their environment and on structure formation as a whole. 
Since their host haloes form from the collapse of rare den- 
sity fiuctuations, miniquasars powered by accreting BHs are 
expected to be strongly clustered and highly biased tracers 
of the underlying dark matter (DM) distribution. We use 
the adaptive mesh refinement (AMR) technique to home 
in, with progressively finer resolution, on the densest parts 
of the 'cosmic web'. Rather than approximating radiation 
fields as isotropic, we study the impact of a point-source of 
X-ray radiation that starts shining in a rare high-cr peak at 
z = 21, before a universal ionizing background is actually 
established. The outline of this paper is as follows. In § 2 
we describe our suite of numerical simulations. § 3 presents 
general results of the simulations, with and without X-rays, 



and on the influence on the thermal and ionization state 
of the IGM of a non-uniform radiation field in the optically 
thin limit. We discuss the distance-dependent radiative feed- 
back effect of a miniquasar on the amount of primordial gas 
that can cool and collapse, thus becoming available for star 
formation, in § 4. Finally, we present our conclusions in §5. 



2 SIMULATIONS 

High-resolution hydrodynamics simulations of early struc- 
ture formation in ACDM cosmologies are a powerful tool to 
track in detail the thermal and ionization history of a clumpy 
IGM and guide studies of early reheating. We have used 
a modified version of enzo, an adaptive mesh refinement 
(AMR), grid-based hybrid (hydro-)- N-body) code developed 
by Bryan & Norman (see http://cosmos.ucsd.edu/enzo/ I to 
solve the cosmological hydrodynamics equations, study the 
cooling and collapse of primordial gas in DM haloes, and 
simulate the large-scale effect of a miniquasar turning on 
at very high redshift. All the results shown below assume a 
ACDM world model with parameters = 0.3, Ha ~ 0.7, 
h = 0.7, Qt = 0.05, as = 0.85, and n = 1. 

The primordial distribution of gas and DM is initialized 
in a comoving box of 1 Mpc on a side using the linear power 
spectrum of Eisenstein and Hu (1999). We first identify in a 
pure N-body simulation the Lagrangian volume of the most 
massive DM halo ai z = 25. At this time its total mass is 
7x10^ Mq , corresponding to a 3.5 sigma peak in the density 
fiuctuation field. Because its mass is comparable to the mass 
threshold for gas cloud formation by molecular cooling, this 
halo will likely harbor the first collapsed baryonic object in 
our simulated volume. It is this halo that we fiag as the 
host of the first massive BH. We then generate new initial 
conditions centered around this density peak, consisting of 
two nested 128^ static grids, of which the inner one covers 
the central 0.5 Mpc volume. The DM density field is also 
sampled with 128^ particles in the inner region, leading to 
a mass resolution of ttidm = 2000 Mq. This ensures that 
haloes above the cosmological Jeans mass are well resolved 
at all redshifts 2;<21. Atz = 15, the five most massive 
haloes in the box have between 5 x 10^ and lO"* DM particles. 

During the evolution from z — 99 to z — 15, refined 
(child) grids are introduced with twice the spatial resolu- 
tion of the coarser (parent) grid. The AMR is restricted 
to the inner 0.5 Mpc and is triggered when a cell reaches 
a DM overdensity (baryonic overdensity) of 2.0 (4.0). This 
low overdensity refinement criterion has been shown to yield 
comparable results in equivalent simulations with enzo and 
the smoothed particle hydrodynamics (SPH) code GAD- 
GET over the entire range of the dark halo mass function 
(O'Shea et al. 2005a). The region of interest is allowed to 
dynamically refine further to a total of 8 levels on a 128'^ 
top grid, resulting in a maximum dynamic range (ratio of 
box size to the smallest spatial scale that can be resolved) 
of 32,768 or a maximum resolution of 30 pc (comoving). At 
the end of a typical simulation, the code uses > 2.2 x 10^ 
computational grid cells on > 9800 grids, with a median 
number of cells per grid of 768. 

The spatial distribution of 128^ x 15/8 coUisionless dark 
matter particles determines at every time-step the large- 
scale gravitational field in the box. We evolve the non- 
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Table 1. Summary of simulations. 



Name Spectrum Energy range Zon Lifetime Description 



NoBH — - — — Reference simulation without miniquasar 

PL power law 0.2-10 keV 21-17.5 45 Myr Fiducial simulation with miniquasar 

PLlong power law 0.2-10 keV 21-15.5 90 Myr Extended miniquasar lifetime 

PLhard power law 0.4-10 keV 21-17.5 45 Myr Harder spectrum 

MCD multicolour disc -|- power law - 21-17.5 45 Myr Includes Lyman- Werner H2-dissociating flux 
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Figure 1. The three types of miniquasar spectra used in this 
paper. 'PL' (solid) and 'PLhard' (dashed) are simple power laws 
with a = I and lower-energy cutoffs of 0.2 and 0.4 keV, respec- 
tively. 'MCD' (dot-dashed) is a combination of multicolour disc 
and a = 1.2 power law. The two components are shown separately 
in light grey. In the MCD case only we have included radiation in 
the Lyman- Werner band (11.2 — 13.6 eV): photons between 13.6 
eV and 0.2 keV are assumed to be absorbed within the emitting 
host halo. All spectra have been normalized to the Eddington lu- 
minosity of a 150 M0 black hole with radiative efficiency e = 0.1. 



equilibrium rate equations for nine species (H, H"*", H~, e, 
He, He^, He"'"'', H2, and H2^) in primordial self-gravitating 
gas, including radiative losses from atomic and molecular 
line cooling, Compton heating by X-rays and Compton cool- 
ing by the cosmic background radiation. The rate coefficients 
for the reaction network for these species are primarily those 
described in Abel et al. (1997): the H2 cooling function is 
that of Lepp & ShuU (1983). We impose the constraint that 
the minimum temperature attainable by radiative cooling is 
that of the cosmic microwave background. 

A zero-metallicity progenitor star of mass m* = 
150 Mq turning on at 2: = 25 in the host halo will emit 
> 10'^* photons above 13.6 eV over a lifetime of 2.3 Myr 
(Schaerer 2002), about 10 such photons per halo baryon. 
The ensuing ionization front will overrun the host halo, pho- 
toevaporating most of the surrounding gas (Whalen, Abel, 
& Norman 2004). Gas accretion on to the BH remnants of 
the first stars that create Hll regions may have to wait for 
new cold material to be made available through the hierar- 
chical merging of many gaseous subunits. The simulations 
show that, hy z — 21, our massive BH has been incorpo- 
rated into a larger DM halo that contains 7 times more gas 



than the original host: it is at this time that we turn on the 
miniquasar radiation field. Our miniquasar is powered by a 
rriBH = 150 Mq hole accreting at the Eddington rate and 
shining for 1 (2 in one of our 7 simulations) Salpeter time- 
scale, ts = 450e/(l — e) Myr, with a radiative efficiency of 
e = 10 per cent. Its exponentially growing mass and lumi- 
nosity are recomputed at every time-step, and its trajectory 
is followed by ffagging the DM particle corresponding to the 
maximum density in the host halo as our BH. The location 
of this flagged particle is determined at every time-step and 
used as the origin of an isotropic 1/r^ radiation fleld: from 
2 = 21 to 2 = 15 our BH typically moves by less than 6 kpc 
(comoving). Note that, because of the considerable mass dif- 
ference between our DM particle resolution (2000 M0) and 
the initial mass of the black hole (150 Mq), we are not able 
to accurately track the true trajectory of the hole. 

The photon energy distribution of putative high-z mini- 
quasars is uncertain. The spectra of 'ultraluminous' X-ray 
sources in nearby galaxies appear to require both a soft com- 
ponent (well flt by a cool multicolour disc blackbody with 
fcTmax — 0.15 keV, which may indicated intermediate-mass 
BHs; Miller & Colbert 2004) and a non-thermal power-law 
component, Le oc E~" , of comparable luminosity and slope 
Q ~ 1. Here we have run simulations with three different 
spectra: one ('PL') assuming a simple power-law miniquasar 
spectrum with a = 1 for photons with energies in the range 
0.2 — 10 keV; one ('PLhard') with a similar spectrum, but 
with a low-energy cutoff of 0.4 keV; and one ('MCD') with 
luminosities equally divided between a multicolour disc com- 
ponent and a power law with a — 1.2, both with a 0.2 
keV low-energy cutoff. In the multicolour disc component, 
each annulus of a thin accretion disc is assumed to radi- 
ate as a blackbody with a radius dependent temperature, 
T(r) oc r~^'^*, and the temperature of the innermost portion 

— 1/4 

of the disc is related to the mass of the BH as Tin cx m^jj 
(e.g. Makishima et al. 2000). All spectral energy distribu- 
tions are normalized to the Eddington luminosity, which 
is exponentially increasing during the active phase of the 
miniquasar. A plot of the different input spectra is shown 
in Fig. For the MCD case only, we included in addition 
to the X-rays an H2-dissociating flux in the Lyman- Werner 
band (LW, 11.2 — 13.6 eV), with an intensity determined by 
the Raleigh- Jeans tail of the multicolour disc. The implicit 
assumption for all of these spectra is that photons with en- 
ergies between 13.6 eV and 0.2 keV are absorbed within the 
miniquasar host halo. For this reason our simulation is not 
strictly valid within the host halo, and we have excised it 
from all of the following analysis. 

In the pre-reionization Universe, when the IGM is pre- 
dominantly neutral, soft X-ray photons will be absorbed as 
they photoionize hydrogen or helium atoms. For a mixture of 
H and He with cosmic abundances, the effective bound-free 
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Figure 2. A 3D volume rendering of the IGM in the inner 0.5 Mpc simulated box at z = 21 and z = 15.5. Only gas with overdensity 
4 < (5 < 10 is shown; the locations of dark matter minihaloes are marked by spheres. The size and grey scale of the spheres indicate halo 
mass. At these epochs, the halo finder algorithm identifies 55 {z = 21) and 262 (z = 15.5) minihaloes within the simulated volume. 



absorption cross-section can be approximated to an accu- 
racy of 30 per cent in the range 50 eV < hv < 2 keV as 
ffbf « 4 X 10"^° cm^ {hv/Q.l keV)"^. The mean free path of 
ionizing radiation in the neutral medium with overdensity 
5 = pI pb is then 



5kpc 



20 



hv 



0.1 keV 



(1) 



Throughout the simulation box, gas at the mean density 
will be transparent to photons with energies > 0.2 keV and 
we take advantage of this by working in the optically thin 
approximation. At 2 = 21 even the most massive haloes 
in our simulation have total hydrogen columns below 10^^ 
cm~'^, i.e. are transparent to photons above 0.7 keV. Since 
a power-law spectrum with v Ii, ^ const is characterized by 
equal power per logarithmic frequency interval, photoelec- 
tric absorption by intervening haloes will not significantly 
attenuate the ionizing energy flux: this is then only a func- 
tion of distance from the miniquasar, and grows linearly with 
the mass of the hole in the PL and PLhard cases. Photo- 
ionizations and photodissociations couple the radiation field 
to the hydrodynamics through a heating term in the en- 
ergy conservation equation and source and sink terms in 
the species abundance equations (Anninos et al. 1997). The 
radiative heating Hj and photoionization and photodissoci- 
ation Tj rate coefficients are given by 



hv — hvn i , 
' ■ dv 



hv 



(2) 
(3) 



where Ojiv) is the cross-section and t'o.j is the frequency 
threshold for the j^^ reaction, and liy) is the intensity of 
the radial radiation fleld. 



3 GENERAL RESULTS 
3.1 'NoBH' simulation 

A reference simulation with no radiation source ('NoBH') 
was run for comparison. The clustered structure around the 
selected high-cr peak of the density field is clearly seen in 
Fig. 121 a 3D volume rendering of the inner 0.5 Mpc box 
at redshifts 21 and 15.5. The flgure shows gas at 4 < 5 < 
10, with the locations of dark matter minihaloes marked 
by spheres colored and sized according to their mass (the 
spheres are only markers, the actual shape of the haloes is 
typically non-spherical). Several interleaving fllaments are 
visible, at the intersections of which minihaloes are typi- 
cally found. We identify the locations of DM minihaloes in 
the simulation volume by employing the HOP halo finder 
developed by Eisenstein & Hut (1998). This algorithm iden- 
tifles haloes by grouping together particles that are asso- 
ciated with the same local density maximum, without en- 
forcing spherical symmetry around the densest particle. We 
consider only haloes containing more than ICQ particles. At 
z — 21, 55 haloes in our volume satisfy this criterion, by 
z = 17.5 this number has grown to 149, and hy z — 15.5 to 
262 haloes. At this epoch, only four haloes have reached the 
critical virial temperature for atomic cooling, Tvir = 10* K, 
where 



1.7 X 10* 



JL 
1.2 

Avir 



107/1-1 M 



2/3 



1 1/3 



QAiiz) 187r2 



20 



K, 



(4) 



p. is the mean molecular weight, and Avir is the density 
contrast at virialization. 

The primordial fractional abundance of H2 in the IGM 
is small, XH2 ~ 2 x 10~®, as at 2; > 100 H2 forma- 
tion is inhibited because the required intermediaries, ei- 
ther Hg""" or H~, are destroyed by cosmic microwave back- 
ground (CMB) photons (e.g. GaUi & Palla 1998). Most of 
the gas in the simulation therefore cools by adiabatic ex- 
pansion. Within collapsing minihaloes, however, gas den- 
sities and temperatures are large enough that H2 forma- 
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Figure 3. Left: Fraction of cold and cold+dense gas within the virial radius of all haloes indentified at z = 17.5 with T^ir > 400, as a 
function of halo mass. Triangles: fc, fraction of halo gas with T < 0.5Tvir and 5 > 1000 that cooled via roto-vibrational transitions of 
H2. Diamonds: /cd, fraction of gas with T < O.STvij and p > 10^^ Mq Mpc~^ that is available for star formation. The straight lines 
represent mean regression analyses of fc and fed with the logarithm of halo mass. Right: Mass-weighted mean H2 fraction as a function 
of virial temperature for all haloes at 2 = 17.5 with Tvir > 400 K and f^d < 0.1 {empty circles) or f^d > 0.1 {filled circles). The straight 
line marks the scaling of the temperature-dependent asymptotic molecular fraction. 



tion is catalyzed by H~ ions through the associative detach- 
ment reaction H+H~ —> H2-|-e~, and the molecular frac- 
tion increases at the rate dxn^/dt oc XeTimT^ii^ , where Xe 
is the number of electrons per hydrogen atom. For Tvir < 
a few thousand kelvins the virialization shock is not ioniz- 
ing, the free electrons left over from recombination are de- 
pleted in the denser regions, and the production of H2 stalls 
at a temperature-dependent asymptotic molecular fraction 
XH2 ~ 10~* TviV^ ln(l-|-t/trec), where Uec is the hydrogen re- 
combination time-scale (Tegmark et al. 1997). A typical H2 
fraction in excess of 200 times the primordial value is there- 
fore produced after the collapse of structures with virial tem- 
peratures of order 10^ K. This is large enough to efficiently 
cool the gas and allow it to collapse within a Hubble time 
unless significant heating occurs during this phase (Abel et 
al. 2000; Yoshida et al. 2003). 

Fig. (left panel) shows the fraction of cold gas within 
the virial radius as a function of halo mass for all the haloes 
identified at redshift 17.5. Note that the presence of halo 
masses below the HOP selection limit of 2 x 10® Mq is 
caused by differences between the extent of a halo as de- 
fined by the HOP algorithm and the virial radius measured 
by the spherically-averaged radial profile analysis. Following 
Machacek et al. (2001), we define fc as the fraction of gas 
with temperature < 0.5 Tvir and density > 1000 times the 
background (this is the halo gas that is able to cool below the 
virial temperature because of H2), and fed as the fraction 
of gas with temperature < 0.5 Tvir and (physical) density 
> lO'^® Mq Mpc~^ (this is the self-gravitating gas available 
for star formation). As in Machacek et al. (2001), we find 
that both fc and fed are correlated with the halo mass, with 
the fraction of cold+dense gas increasing less rapidly than 
fc with Mhaio. The mass threshold for significant baryonic 
condensation (non-zero fed) is approximately 5 x 10® M0 at 



these redshifts (Haiman, Thoul, & Loeb 1996). Also depicted 
in Fig. 13 (right panel) is the mass- weighted mean molecular 
fraction of all haloes with Tvir > 400 K. Filled circles repre- 
sent haloes with fed > 0.1, while open circles represent the 
others. The straight line marks the scaling of the asymp- 
totic molecular fraction in the electron-depletion transition 
regime. Haloes with significant cold-|-dense gas fraction have 
larger (xhs) than those shown in fig. 3 of Yoshida et al. 
(2003). This is due to the high resolution attainable by an 
AMR code compared to SPH. The maximum gas density 
reached at redshift 15 in the most refined region of our sim- 
ulation is 4 X 10® cm~^ (corresponding to an overdensity of 
3 X 10*): within this cold pocket the excited states of H2 
are in LTE and the cooling time is nearly independent of 
density (e.g. Lepp & ShuU 1983). 

At z — 17.5, ~ 90 per cent of the simulated volume 
contains gas with T < 10 K, as shown by the differen- 
tial temperature distributions in Fig. |1| The distribution of 
the entropy parameter k = kTn~^^^ is even more strongly 
peaked than the temperature distribution, as entropy is con- 
served during infall on to non-linear structures or Hubble 
expansion, which heat/cool the gas adiabatically. After de- 
coupling from the CMB, the IGM temperatures drops as 
T{z) ~ 2.73(1 + ^d)[(l + 2)/(l + ^d)]'K for z < Zi ^ 150, 
giving rise to an entropy ffoor of 0.04 eV cm'^, independent 
of redshift. Only gas that is shock heated during accretion 
on to DM haloes or that cools radiatively via H2 can depart 
from this peak. Shock heating produces the extended tail 
out to K = 10 eV cm^ at a volume fraction in the range 
fv = 10"® - 10~*. At these epochs, we find a collapsed gas 
mass fraction of /a/ = 2.5 per cent. The fraction of baryons 
that gets shock heated (k > 0.1 eV cm^) is /ji/ = 6.7 per 
cent. Radiative cooling processes affects such a tiny fraction 
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Figure 4. Distributions of temperature T (binsize AlogT = 0.04), 'entropy' k (AlogK = 0.08), electron fraction (AlogXe = 0.07), 
and molecular liydrogen fraction x^i^ (^log xjjj = 0.10) (clockwise from upper left). Dotted lines: z = 17.5 NoBH run without miniquasar 
feedback. Dark solid lines: z = 17.5 fiducial simulation, in which the miniquasar shines for one Salpeter time (2 = 21 — 17.5). Light solid 
lines: z = 15.5 PLlong run, in which the miniquasar shines for two Salpeter times (z = 21 — 15.5). Dashed lines: z = 17.5 PLhard run, 
in which the low-energy cutoff of the miniquasar spectrum has been increased from 0.2 to 0.4 keV. The volume- averaged temperature, 
entropy, electron fraction, and molecular fraction are (8K,0.04eV cm^,lA X 10"'', 1.9 X 10"*^), (2800K,17eV cm^, 0.03, 3.9 X lO"'^), 
(6150K,49eV cm^, 0.07, 9.3 x lO'^), and (383K,2.2eV cm^, 0.005, 6.3 x lO'^), respectively 



of the volume that it remains outside the plotted range. 
Due to the lack of ionizing sources the medium is almost 
completely neutral, with a distribution of electron fractions 
sharply peaked at — 1.5 x 10~*. Note that the distri- 
bution of molecular hydrogen fraction, xn^, is considerably 
broader than that of Xe- 

Another quantity of interest is the gas clumping factor 
C — (n^)/(n)^, where the brackets denote volume- averaged 
quantities. The shorter radiative recombination time-scale, 
tree ~ {a{ne)C)^^ of a clumpy medium relative to a ho- 
mogeneous one increases the total number of photons per 
baryon required for reionization. We have calculated the av- 
erage clumping factor within our simulation box of volume 
as 

C = CiGM + Ch = -; TT7T + 1 T^TTT) (5) 

<S<70 ^ ' i>70 ^ ' 



where Ui is the number density of hydrogen in the i"" grid 
cell of volume Vi, (p-h) is the background hydrogen den- 
sity, and the two sums are taken over all grid cells with 
baryonic overdensity smaller and larger than 70. We use 
this density threshold - corresponding to the mean baryonic 
overdensity at the virial radius of our haloes - to differen- 
tiate dense gas belonging to virialized haloes (Ch) and dif- 
fuse uncoUapsed 'intergalactic' material (Cigm). Gas within 
haloes only contributes to recombinations when it is pho- 
toionized: self-shielding of a cold gaseous disc and photoe- 
vaporation of halo gas by an external radiation field all make 
the contribution of dense material with 5 > 70 to the over- 
all clumping very uncertain (e.g. Haiman et al. 2001). In 
Fig. |S]we have plotted the evolution with redshift of Cigm 
and Ch in the NoBH simulation. The clumping factor in- 
creases rapidly with time as structure formation progresses, 
with Cigm ~ 2.65 exp[0.2(21 — z)] over the range in redshift 
we probe. 



^ The effect of radiative cooling can be observed in the 'phase 
diagram' of Fig.EJas the bluish tail at <5 > 100, T < 1000 K, and 
volume fraction fy < 10~*. 
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Figure 5. Clumping factor {{n^) / {n)'^) as a function of redshift for the NoBH, PL, and PLlong simulations. Left: Evolution of the 
clumping factor associated with diffuse 'intergalactic' (5 < 70) material. Right: halo (5 > 70) clumping factor. Note that this is not 
reduced in the presence of X-ray heating. 
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Figure 6. Two-dimensional distribution of gas temperature versus baryonic overdensity at z = 17.5. The color coding indicates the 
fraction of the simulated volume at a given (5, T). Left: NoBH run. Right: PL run. 



3.2 Heating and ionization by the miniquasar 

X-ray radiation from the miniquasar partially ionizes most 
of tlie gas in tiie simulated volume botii by direct photoion- 
ization of H I and He I and indirectly by coUisional ionization 
of H I by tlie fast photolectrons. We liave adopted the fitting 
formulae of Shull & van Steenberg (1985), cast as a function 
of the hydrogen ionized fraction, x = njj+/nH, to deter- 
mine the fraction of energy of a photoelectron deposited as 
heat versus further ionization. For photon energies above 0.2 
keV, the ratio of the H I /He I photoionization cross-sections 
drops below 4 per cent; since the primordial ratio of helium 
to hydrogen is about 8 per cent, the photoionization and 
heating rates are dominated by helium absorption, which 
exceeds the hydrogen contribution by ~ 2 : 1. While He I is 
the main source of hot primary photoelectrons, however, it 
is H I that undergoes the bulk of secondary ionizations. A 
primary nonthermal photoelectron of energy E — 1 keV in a 



medium with residual ionization fraction (from the recom- 
bination epoch) X = 2 X 10~* will create over two dozens 
secondary electrons, depositing a fraction /i — 0.37 of its 
initial energy as secondary ionizations of Hi , /2 = 0.05 
as secondary ionizations of He: , and /a = 0.13 as heat. 
The time-scale for electron-electron encounters resulting in 
a fractional energy loss / = AE/E, 



to 



140 yr Ef 



l + z 
20 



In A 
"20" 



(6) 



(where E is measured in keV), is typically much shorter 
than the electron Compton cooling time-scale: the primary 
photoelectron will therefore ionize and heat the surround- 
ing medium before it is cooled by the CMB. The fraction 
of primary energy going into heat increases gradually with 
X, reaching /a = 0.5 at a: = 0.04. Secondary ionizations en- 
hance the rates for ionization of Hi and He I (but not for 



© 2005 RAS, MNRAS OOO.ITIITsl 



8 Kuhlen & Madau 



Hen , ShuU & Van Steenberg 1985), and introduce an addi- 
tional coupling between the abundance rate equations. The 
rate coefficients for Hi and He I ionizations can be written 
as 

rH: + /.(.)(/fH: + ^i/H.)^ (7) 
and 

rH. + /.(.)(i^H. + ^/fH:)^, (8) 

respectively. We self-consistently evolve the hydrogen ion- 
ization fraction x and the hydrogen and helium abundances 
throughout the simulations (note that the ShuU & van Steen- 
berg fitting formulae were derived assuming equal ionization 
fractions of hydrogen and singly-ionized helium) . 

The two-dimensional distribution of gas overdensity and 
temperature at z = 17.5 is shown in Fig.Elfor the NoBH and 
PL simulations. The color coding in this phase diagram indi- 
cates the fraction of the simulated volume at a given (5, T). 
Most of the gas in the NoBH run lies on the yellow line 
(k = KTn~^^^T = const) representing the initial adiabat. 
At low overdensities the temperature either drops because of 
Hubble expansion or rises because of adiabatic compression 
until the gas is shock heated (red and green swath) to virial 
values, T = 10'^ — 10* K. At higher densities, the blue cooling 
branch follows the evolutionary tracks in the temperature- 
density plane for spherically collapsing clouds (Yoshida et 
al. 2003). H2 line emission lowers the temperature down to 
~ 100 K, the minimum value attainable by molecular cool- 
ing. The onset of the gravitational instability can further 
compress the gas and cause a modest rise in temperature 
again (Abel et al. 2000). 

After shining for a Salpeter time-scale, the miniquasar 
has heated up the box to a volume- averaged temperature of 
2800 K. The mean electron fraction and entropy are now 0.03 
and 17 eV cm^ (see Fig.|3|l: farther than 20 kpc (comoving) 
from the source hydrogen is never ionized to more than 30 
per cent. Gas near the miniquasar is heated above 10* K 
and quickly cools down via efficient atomic processes: the 
green finger at log 5 = O,log(r/K) > 4 represents baryonic 
material in this phase. The increased electron fraction and 
gas temperature boost the gas-phase H2 production, which 
occurs on a time-scale, 

that is much shorter than the hydrogen recombination time, 

W = 9 X 10^ Myr (^^) T^^S'' (^i±i) (10) 

(here T3 = T/IO^ K). The volume -averaged molecular frac- 
tion is now 20 times larger than the primordial value, with 
denser filaments in the IGM (5 ~ 10 — 20) being traced out 
by a smaller electron fraction and by a molecular fraction 
in the range 10~* — lO"'^ (see Fig. |7|l. This large molecular 
fraction will hinder the buildup of a uniform UV photodis- 
sociating background, as the maximum optical depth of the 
IGM in the H2 LW bands can now exceed 10 (cf. Ricotti, 
Gnedin, & ShuU 2001; Haiman et al. 2000). It wiU also pro- 
mote H2 radiative cooling in filaments on a time-scale. 




10° 10' lO"" 10^ lO"* 10^ 
6 



Figure 8. Mean H2 cooling time (in units of H~^) versus gas 
overdensity, for the NoBH (solid line) and PL (dashed line) runs 
at z = 17.5. The thin lines show departures of one root mean 
square (((t^^^i) — (tcool)'^)^^'^) from the mean. Gas below the thin 
horizontal line can cool down in a Hubble time. 

« 650Myr (^H^) x^^S'^ (t)'' ^^'^ 

(Machacek et al. 2001) that can be shorter than the Hubble 
expansion time, 1/H ~ 285 Myr [(1 -I- 2)/20]"^/^, and the 
Compton cooling time for mostly neutral primordial gas, 
tc = SMyrx-^ [(1 + 2)/20]"*. In Fig. Hwe plot the mean 
H2 cooling time in units of the expansion time-scale 
versus overdensity, sd, z — 17.5. In the NoBH simulation only 
gas with 5 > 600 can cool in a Hubble time. This critical 
cooling density is lowered to 5 ~ 40 in the PL simulation, 
due to the enhanced H2 fraction. Molecular cooling can then 
affect the dynamics of baryonic material before it has fallen 
into the potential well of DM haloes and virialized. Gas in 
the simulation box is now able to collapse and cool more 
rapidly within subsequent star-forming minihaloes, without 
the need to produce much additional H2. 

The global environmental impact of the miniquasar is 
illustrated in Fig. |5] where we plot in the top panel the vol- 
ume (fv) and mass fraction (/m) as a function of overden- 
sity in the NoBH and PL runs. The bottom panel shows the 
corresponding percentage change from the NoBH to the PL 
case (100 X (/pl — /noBh)//noBh). Three prominent features 
are clearly seen in this figure: 1) over most of the simulation 
volume, the net effect of X-rays is to reduce gas clumping 
due to the smoothing of sheets and filaments in the 'cos- 
mic web' by gas pressure. This 'Jeans smoothing' removes 
gas from intermediate S's and puffs it up to low overdensi- 
ties (A/ > for 1 < 5 < 10), thereby filling in underdense 
regions (A/ < for S < 1);^ 2) the suppression of bary- 
onic infall and the photoevaporation back into the IGM of 
some of the gas already incorporated into haloes both lower 

^ In the miniquasar simulations the IGM clumping factor remains 
approximately constant, Cigm — 3, from z = 21 to 2 = 15 (see 
the left panel in Fig. |5J. Because of the decreasing background 
density with time, the IGM recombination time-scale at a given 
temperature therefore increases as (1 -|- z)~^ over this redshift 
interval. 
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Figure 7. From top to bottom: Projected temperature, electron fraction, and molecular fraction distributions in the simulation box for 
runs NoBH (left) and PL (right) at z = 17.5. The bottom panel shows the projected gas density in a smaller 127 comoving kpc region 
(dotted outline in top panels) to highlight the effect of Jeans smoothing. 
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the gas mass fraction at 20 < 5 < 2000; and 3) enhanced 
molecular cooling increases the amount of dense material 
at (5 > 2000. While the vast majority of the baryons, both 
in volume and mass, are heated up, the densest gas at the 
centres of haloes is actually cooled, due to the X-ray cataly- 
sis of molecular hydrogen. The feedback effect of X-rays on 
the formation and cooling of pregalactic minihaloes will be 
discussed in more details in the next section. 

In the PLlong simulation we let the miniquasar shine 
for two Salpeter times, from 2: = 21 to 2 = 15.5, leading 
to a final black hole mass of 1100 Mq. The increased in- 
tensity and duration of the radiation further enhances the 
effects described above (see Fig.^J. At 2 = 15.5, the volume 
averages in the box are: (T) = 6150 K, (k) — 49 eVcm^, 
(xe) = 0.07, and (xHa) = 9.3 x 10"^. The additional heat- 
ing, however, does not further decrease the IGM clumping 
factor, and the halo clumping factor is still not affected (see 
FigEl • The increase in {x-n^ ) by more than a factor of two 
is the most significant difi'erence between the PL to PLlong 
simulation. This will make it even harder for subsequent 
generation of sources to establish a uniform Lyman-Werner 
background. Note that if early black holes undergo a period 
of super-Eddington mass growth (Volonteri & Rees 2005; 
Haiman 2004), a possible subsequent luminous miniquasar 
phase, powered by a more massive hole, might have feed- 
back effects more similar to, or even exceeding those of our 
PLlong simulation. 

In order to probe the dependence of our results on the 
hardness of the radiation field, we ran the PLhard simula- 
tion, in which the lower energy cutoff was raised to 1^111111=0.4 
keV. Normalizing the total energy output to LEdd leads to 
a specific radiation intensity which is 20 per cent higher (for 
Q = 1) compared to the PL simulation. The heating rate, 
on the other hand, scales roughly as v^f^. A factor of two 
increase in t-min thus leads to a factor of eight decrease in 
the heating rate, and we expect the feedback from the mini- 
quasar to be much weaker. This expectation is confirmed by 
the PLhard simulation, as shown by the dashed histogram in 
Fig- ID At 2 = 17.5 the volume-averaged temperature is only 
380 K, with (xe) = 0.005 and (xHa) = 6.2 x 10"^. By the 
same arguments we would expect the miniquasar feedback 
to strengthen if the spectrum extended to lower energies, 
although the optically thin assumption would then break 
down. 



4 X-RAY RADIATIVE FEEDBACK: POSITIVE 
OR NEGATIVE? 

The impact of the radiation produced by Population III 
stars and their remnants on H2 chemistry and the proper- 
ties of star-forming gas in low-mass haloes has been studied 
by many authors. Molecular hydrogen is fragile and easily 
photodissociated by soft UV photons below 13.6 eV (e.g. 
Haiman, Rees, & Loeb 1997; Ciardi, Ferrara, & Abel 2000; 
Machacek et al. 2001; Glover & Brand 2001). While this 
negative feedback may suppress gas cooling and collapse, 
the formation of H2 molecules is catalyzed by free electrons, 
and any process that increased their abundance would boost 
the abundance of H2 as well. Fossil Hll regions (Ricotti et 
al. 2001) and X-ray photons may provide such positive feed- 
back (Haiman et al. 2000; Venkatesan, Giroux, & Shull 2001; 




Figure 9. Top: volume (black lines) and mass {gray lines) frac- 
tions of gas overdensity 5 = p/pb, in the NoBH (solid lines) and 
PL (dashed lines) simulations at z = 17.5 (binsize Alog(5 = 0.09). 
Bottom: Percentage change, 100 X [/(PL) - /(NoBH)]//(NoBH)] 
in gas volume/mass fraction as a function of S. 



Machacek et al. 2003; Glover & Brand 2003; Cen 2003). The 
issue of whether positive feedback processes dominate over 
negative feedback is still open. All we offer here is an as- 
sessment of the effect of X-rays in the absence of a strong 
photodissociating fiux. This condition would naturally be 
fulfilled if accreting intermediate-mass BHs were the end- 
product of the first episodes of Population III star formation, 
as miniquasars are much more efficient sources of radiation 
than their stellar-progenitors (Madau et al. 2004). As al- 
ready mentioned in § 2, our MCD run includes in addition 
to the X-rays a LW H2-dissociating fiux with an intensity 
determined by the Raleigh- Jeans tail of the multicolour disc, 
Ilw{z = 21) = 10~^^ ergs cm~^ s~^ Hz~^ sr~^ at a comov- 
ing distance of 7.2 kpc, which is 10 (3) times lower than the 
0.2 keV (1 keV) X-ray fiux (cL Machacek et al. 2003). We 
find the results of the MCD simulation to be very similar to 
those of the PL run, and it is the latter that we analyze in 
detail below. 

Fig. 1101 depicts the change in halo gas content between 
the PL and NoBH simulations as a function of halo mass 
at two different epochs, 2 = 17.5 (when the miniquasar 
stops shining) and 2 = 15.5: haloes in the two simulations 
have been matched by location and DM mass. Different sym- 
bols represent total, cold, and cold-l-dense gas mass both for 
haloes closer than 75 comoving kpc from the radiation source 
and for those farther away. X-ray heating lowers the total 
gas mass in nearly all virialized DM clumps, both by pho- 
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Figure 10. Variation in halo gas mass between the PL and NoBH simulations as a function of halo total mass (as measured in the NoBH 
run) at two different epochs. Symbols show log[Mi(PL)/Mi(NoBH)], with i = g,c,cd for total gas, cold gas, and cold+dense gas mass. 
Filled symbols represent haloes within 75 comoving kpc of the miniquasars, empty symbols the others. Haloes with |AlogA/(. cdl > 2 are 
shown as arrows at the top and bottom of the plot. 



toevaporating baryons that were already incorporated into 
haloes and by suppressing baryonic infall into newly form- 
ing structures. This decrease in gas mass is more apparent 
below Mhaio ~ 10® M0, and it is the strongest in pregalac- 
tic clouds close to the miniquasars that have been exposed 
to the highest level of X-ray flux. Only a small fraction of 
the gas accreted by haloes in the NoBH run is retained by 
the host gravitational potential in the PL run. One of the 
potential wells that experiences the largest decrease in gas 
content is fairly massive, Mhaio = 10® M©, is located only 12 
kpc (comoving) away from the miniquasar, and has a total 
gas fraction of 0.7 per cent, more than 20 times lower than 
in the NoBH simulation! 

Population HI stars, however, can only form from gas 
that is able to cool and condense. Fig. 1111 shows the mass 
weighted mean H2 fraction inside haloes in the PL simula- 
tion at z — 17.5 (cf. right panel of Fig. Overall (a;H2) 
increases by slightly more than one order of magnitude, in- 
dependently of Tvir- Compared with the NoBH case, haloes 
with fed > 0.1 can be found at slightly lower T^ir- None of 
the haloes with Tvir < 800K has any cold-|-dense gas, even 
though their molecular fraction is quite high, {xti2) ^ 10"'^. 
These haloes are the least massive haloes in our simula- 
tion and were formed only after the miniquasar began to 
shine. Continuous X-ray heating has prevented the gas from 
cooling and falling into the halo. As a consequence their 
mean gas mass fraction (Mgas/Mtot) is only 0.08 and their 
mean central temperature exceeds their virial temperature 
by a factor of 3. Haloes with larger Tvir, in contrast, formed 
earlier and their gas has contracted to densities where H2 
cooling is able to overcome X-ray heating. 

Evidence for X-ray enhanced cooling is seen in many 
haloes above 2 x 10^ M0. The largest relative increase 
in the amount of cold material occurs in the mass range 
2 X 10''' < A/halo < 10® M0, where the boosting eflect can 
exceed 1-2 orders of magnitude! In this sense. X-ray pre- 



Table 2. Summary of miniquasar feedback on halo total, cold, 
and cold+dense gas mass, at z = 17.5. All numbers are per cent 
change from NoBH to PL case. 

O.SMpc box < 75kpc > 75kpc > ISOkpc 

total -20.5 -57.4 -10.5 -8.2 

cold -1-19.4 -43.4 -f43.5 -f49.3 

cold+dense +5.7 -52.1 +29.3 +38.4 

heating appears to somewhat decrease the threshold mass 
required for efficient gas cooling and star formation, albeit 
not by a large factor. In more massive peaks most baryons 
are already cold in the absence of X-rays (see Fig. |2J , and 
positive feedback can only promote a little additional cool- 
ing. Yet in many haloes in the proximity of our miniquasar 
the H2boosting eflect of X-rays is too weak to overcome 
heating, and the cold and dense gas mass actually decreases. 
Table |5| summarizes the miniquasar's feedback effect on the 
total, cold, and cold+dense gas mass. We find evidence for 
a global positive feedback in minihaloes more than 75 co- 
moving kpc away from the miniquasar, as the total mass in 
cold and cold+dense baryons increase by 30 and 44 per cent, 
respectively. Farther than 150 kpc this positive feedback is 
slightly larger yet, reaching +50 per cent for cold gas. Within 
75 kpc of the miniquasar, however, the feedback is negative, 
as the total mass in cold+dense baryons decreases by more 
than 50 per cent. The total amount of gas mass in haloes is 
reduced at all distances within our simulated volume, with 
haloes closer than 75 kpc losing on average more than half 
their gas mass. These effects are long-lasting, and persist at 
z — 15.5, 40 Myr after the miniquasar has stopped shining 
(see right panel of Fig. IIOII . Reducing the total amount of 
gas in the outer regions of haloes, while raising the amount 
of cold and dense gas available for star formation in the 
centres, could make it easier for these haloes to chemically 
enrich their surroundings. 
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Figure 11. Same as Fig.l^frio/it panel), but for the PL simula- 
tion. The dotted line is the T^-^ line from the NoBH case. 



Overall, however, positive and negative feedbacks nearly 
compensate each other, and the net impact of X-rays on 
the total amount of cold material available for star forma- 
tion within the simulated volume is remarkably mild. By 
z = 17.5, exposure to X-ray radiation has boosted the total 
mass fractions of cold and cold-l-dense gas by only 14 and 6 
per cent, respectively. There is hardly an X-ray 'sterilizing 
effect': star formation can proceed in minihaloes above the 
threshold mass, including those that are near the radiation 
source. In the most massive peaks at 2: = 15.5 the frac- 
tion of gas available for star formation in the NoBH control 
simulation is comparable to the value found when X-rays 
are present. Qualitatively similar conclusions were reached 
by Machacek et al. (2003), who showed that an early X-ray 
background cannot overcome the negative feedback from H2 
photodissociation by the soft UV radiation spectrum of the 
first stellar sources. Note that, at a distance of 50 comoving 
kpc, the X-ray flux from our miniquasars is comparable to 
the strongest X-ray background adopted by Machacek et al. 
(their ex = 10 case). Here we find that, even in the absence 
of a strong LW flux, the radiative feedback from X-rays is 
subtle. It enhances gas cooling in lower-cr peaks that are far 
away from the initial site of star formation, thus decreas- 
ing the clustering bias of the early pregalactic population, 
but does not dramatically reverse or promote the collapse 
of pregalactic clouds as a whole. 

It has been suggested recently by Oh & Haiman (2003) 
that an early X-ray background would establish an entropy 
floor over the entire IGM, thus preventing gas contraction, 
H2 formation, cooling, and the build up of dense cores in 
minihaloes. The implication is a large reduction in the col- 
lapsed gas fraction and a pause in the cosmic star formation 
history, before more massive haloes with Tvir > 10* K (which 
can undergo atomic cooling) start forming. We find little 
evidence for preheating suppressing H2 formation. None of 
our minihaloes exhibits an entropy floor and, as we have 
shown above, the amount of cold and dense gas in the cen- 
tres is actually enhanced in haloes sufficiently removed from 



the miniquasar. Oh & Haiman considered the evolution of 
a pre-heated, isolated, uniform density gas parcel with ini- 
tially primordial H2 fraction, and showed that the entropy 
floor prevented subsequent H2 formation. In our simulations, 
however, the heating source is also a catalyst of H2. This 
causes the H2 cooling time to become shorter than the Hub- 
ble expansion time at 5 > 40 (see Fig. |SJ, and thus gas 
entropy is no longer a conserved quantity inside and in the 
vicinity of minihaloes. 

Fig. 1121 shows spherically averaged mass- weighted pro- 
files of S, T, K, Xe and at 2 = 17.5 of a typical low mass 
minihalo in our simulations. It has a dark matter mass of 
3 X lO'^ M0 (Tvir « 925 K), a virial radius of 127 pc (proper), 
and is located 200 comoving kpc from the miniquasar. In the 
NoBH simulation it has a total gas mass of 5.7 X 10* Mq, 
a cold gas mass of 1.1 x 10* Mq, and no cold-|-dense gas. 
The central regions of this halo have only triggered refine- 
ment down to level 7, one shy of the maximum refinement 
level. It has a central core, where the density profile lev- 
els off at (5 = 7000. The gas temperature rises from about 
100 K to Tvir at the virial radius, and then drops down 
to the central value of 350 K. This drop is possible due 
to the slightly elevated central H2 fraction of 2.7 x 10~*, 
which, however, is not high enough to allow the gas to con- 
dense to the cold-|-dense threshold density of 330 cm""^ in a 
Hubble time. The X-ray fiux from the miniquasar changes 
this picture dramatically. The H2 catalysis has led to an en- 
hancement of almost 2 orders of magnitude in outside 
the virial radius, from 10~^ to 8 x 10~*. Towards the cen- 
tre the increase in temperature and density allows the H2 
abundance to grow further, peaking at a central value of 
XH2 = 6.6 X 10~^. This increase in XH2 causes the gas to 
cool down to 140 K and greatly increases the central over- 
density (5 = 3 X 10^), triggering refinement all the way to 
the maximum refinement level. In the PL simulation the 
external medium is at a temperature of > 3000 K, which 
is much higher than Tvir. Thus no virial shock developes, 
and H2 cooling causes the temperature profile to decreases 
monotonically towards the centre, dropping below Tvir at 
~ 40 pc, a third of the virial radius. Gas that is hotter than 
Tvir cannot accrete on to the halo, and as a consequence 
the outer halo regions show a reduced gas density. The total 
gas mass decreases by 30 per cent down to 4.0 x 10* M©. 
The increased cooling, however, has boosted the amount of 
cold gas by 60 per cent, up to 1.8 x 10* Mq, and even al- 
lowed 6.5 X 10'^ M0 to reach the cold-|-dense threshold. The 
entropy profile reflects the additional cooling in the centre: 
it also decreases monotonically and reaches a minimum of 
2.2 X 10"* eV cm^ , a factor of 30 below the central value in 
the NoBH case (6.4 x 10"^ eV cm^). 

The spatial resolution of our AMR simulations is deter- 
mined by the maximum refinement level, which is usually set 
by computational costs. Mesh cells that have refined down 
to the maximum refinement level can pose a problem. The 
physics of cooling and collapsing gas is inherently unstable: 
Jeans unstable cells should refine further, but the simulation 
does not allow it. This leads to an error in the numerical 
solution of the differential equations governing the system. 
One way to prevent the gas from continuing to trigger re- 
finment is to introduce Artificial Pressure Support (APS) 
to stabilize the density peaks that have reached maximal 
refinement against further collapse. With APS the internal 
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Figure 12. Spherically averaged mass-weighted radial profiles of <5, T, re, Xe, and xn^ of a representative halo a.t z = 17.5. In the left 
panel we compare the NoBH with the PL case, in the right panel the PL with the PLIO case. This halo has a DM mass of 3 X 10^ 
and is located 200 comoving kpc from the miniquasar. See text for discussion. 



gas energy of every mesh cell that reaches maximum refine- 
ment and is Jeans unstable is raised to 10 times the value 
that would stabilize it against gravitational collapse. This 
provides support against further refinement, but comes at 
the cost of introducing an artificially high central pressure 
and temperature inside haloes that have reached maximum 
refinement. Neither applying APS nor letting the simulation 
progress without it is correctly modelling the physics of the 
problem. The only solution is to continue to refine further 
(Abel et al. 2002) or to introduce sink particles (Krumholz, 
McKee, & Klein 2004). In our fiducial simulations we have 
50,000-60,000 mesh cells at the maximum refinement level 
of 8. In order to determine whether the associated numerical 
error affects the conclusions of our study, we have re-run the 
PL simulation from z — 21 to z — 17.5, with a maximum 
refinement level of 10 (for a maximum dynamic range of 
131,072) and APS. In the right panel of Fig.IElwe plot for 
the PL and new 'PLIO' runs a comparison of the five profiles 
described above. As expected the two profiles agree excel- 
lently at radii outside of the central core (> 10 pc). Even 
inside the core, the temperature and profiles agree very 
well. The gas in the PLIO halo, however, has been able to 
reach an overdensity a factor of ~ 3.5 higher than in the 
PL case, and has ~ 50 per cent more cold-|-dense gas. Al- 



though we show here only profiles for one halo, we see similar 
trends for all haloes that reach maximum refinement, and 
so we conclude that the amount of cold-f dense gas in haloes 
in the PL simulation should be viewed as a lower limit. The 
numerical error due to evolving the simulation with cells at 
the maximum refinement level is unlikely to affect the global 
conclusions of this paper. 



5 DISCUSSION 

Active galactic nuclei powered by supermassive holes keep 
the Universe ionized at z < 4, structure the IGM, and proba- 
bly regulate star formation in their host galaxies. Their seeds 
were likely planted at very early epochs through the collapse 
of massive stars. Intermediate-mass holes accreting gas from 
the surrounding medium may shine as miniquasars at red- 
shifts as high as z ~ 20. In this paper we have carried out 
AMR cosmological simulations using enzo to address the 
thermodynamic effect of miniquasars on the IGM at early 
times. 

X-ray radiation from the miniquasar efficiently heats 
the gas in the simulated box to a volume-averaged tempera- 
ture of 2800 K after one Salpeter time {z = 17.5), and 6150 
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K after two {z = 15.5). The main effect of this sharp increase 
in temperature is a reduction of gas clumping in the IGM 
by as much as a factor of 3, due to the Jeans smoothing of 
sheets and filaments in the 'cosmic web' by increased gas 
pressure. This smoothing will lower the number of hydro- 
gen recombinations and thus the number of UV photons per 
baryon required to reionize the universe. Since X-rays are 
more efficient at heating than at ionizing, the free electron 
fraction is never raised to more than ~ 10 per cent, reaching 
a volume average of 3 per cent at 2 = 17.5. Provided that 
either collisions or Lya radiation can couple the spin tem- 
perature to the kinetic temperature of the largely neutral 
gas, it may be possible to observe the X-ray heated region 
in redshifted 21-cm line emission against the CMB (Madau, 
Meiksin, & Rees 1997) with future facilities like the LOw 
Frequency ARray (LOFAR). 

The elevated free electron fraction leads to a strong en- 
hancement of molecular hydrogen, both inside haloes and in 
the more tenuous filaments. The volume- averaged H2 frac- 
tion is raised to = 4 x lO^'', 20 times larger than the 
primordial value. This will delay the buildup of a uniform 
UV photodissociating background by subsequent sources of 
Lyman- Werner photons. We plan to carry out a detailed cal- 
culation of this suppression and its consequences in future 
work. The increased H2 abundance allows gas with over- 
densities (5 > 40 to cool within a Hubble time. As a conse- 
quence the evolution of gas within haloes and filaments is not 
adiabatic, and no entropy fioor is established. While global 
heating suppresses baryonic infall and lowers the gas mass 
fraction at overdensities 5 in the range 20-2000, enhanced 
molecular cooling increases the amount of dense material 
at S > 2000. Thus, while the vast majority of the baryons, 
both in volume and mass, is heated up, the densest gas at 
the centres of haloes is actually cooled down. The largest rel- 
ative increase in the amount of cold material occurs in the 
mass range 2 x 10^ < Mhaio < 10^ M0, where the boost- 
ing effect can exceed 1-2 orders of magnitude. Yet for many 
haloes in the proximity of our miniquasar (< 75 kpc) the 
H2-boosting effect is too weak to overcome heating, and the 
cold and dense gas mass actually decreases. Overall, the ra- 
diative feedback from X-rays enhances gas cooling in lower-cr 
peaks that are far away from the initial site of star formation, 
thus decreasing the clustering bias of the early pregalactic 
population, but does not appear to dramatically reverse or 
promote the collapse of pregalactic clouds as a whole. 

While we have included the essential physics necessary 
to model the radiative feedback from the first miniquasars, 
our study, like any, has its limitations. In the following we 
briefiy address some of these: 

1) We have not considered the environmental impact of 
the progenitor of the black hole. In the case of direct collapse 
of a primordial gas cloud to an intermediate-mass black hole 
the issue of a progenitor star would not arise. The more typ- 
ical formation mechanism of the black hole, however, may 
be the collapse of a Population III star. In this case the copi- 
ous UV radiation emitted during its main-sequence lifetime 
could have significantly altered the initial conditions for our 
simulation. Recent three-dimensional simulations (O'Shea 
et al. 2005b) seem to indicate that the heating and ionizing 
effect from the Population III progenitor are rather short 
lived. In their simulation the gas quickly recombines and 



cools, H2 reforms, and gas is able to collapse to form a sec- 
ond generation star in a neighboring halo only 265 pc from 
the Population III host halo. 

2) The strong increase in H2 abundance catalyzed by the 
X-rays is dependent on the absence of a Lyman-Werner 
H2-dissociating background. As shown by Machacek et al. 
(2003) the positive feedback effects from an X-ray back- 
ground in the presence of a LW background are rather mild, 
and they did not find a global increase in H2 abundance by a 
factor of 20. For this reason our work should be considered 
as a simulation of the first miniquasar in its volume, i.e. 
before other sources establish a LW background. How soon 
after the miniquasar a LW background can be established, 
given the large increase in H2, remains to be seen. 

3) Recently there has been a renewed interest in the 
importance of deuterium hydride (HD) to the cooling of 
primordial gas clouds (Johnson & Bromm 2005; Lipovka, 
Nunez-Lopez, & Avila-Reese 2005; Nagakura & Omukai 
2005). HD becomes more efficient than H2 cooling at T < 
200K, and, given a sufficient HD abundance, can cool gas 
down to below 100 K. In our simulations deuterium chem- 
istry has been neglected. An increased free-electron fraction, 
as produced by our miniquasar, would lead to a catalysis of 
HD molecules, so it is quite possible that HD cooling may 
become important in the cores of our minihaloes. 

4) As the emission spectrum of a miniquasar is uncertain, 
so are the photoheating and photoionization rates. We have 
mainly considered a pure power-law spectrum, ranging from 
0.2 to 10 keV, with a slope of —1. The power-law component 
in our MCD simulation is slightly softer {a — —1.2), but the 
differences between PL and MCD are negligible. 

Note also that we have neglected UV radiation with 
13.6 eV < hv < 200 eV. If a significant amount of the emit- 
ted Lyman-continuum fiux escapes the host halo, it will be 
necessary to follow the radiative transfer in order to fully 
assess the environmental infiuence of the first miniquasars. 

5) We have not taken into account the possibility of a 
feedback effect from the miniquasar's radiation on the black 
hole's mass accretion rate. Instead of assuming a constant 
exponential growth rate, it would be better to determine M 
as a function of time from the amount of cold gas available 
for accretion. Unfortunately our simulations don't resolve 
the length scale of the accretion ffow, and a self-consistent 
treatment of the mass accretion rate will have to wait for 
higher resolution studies. 
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